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Abstract 

We interpret within the phenomenological MSSM (pMSSM) the resuhs of SUSY 
searches published by the CMS collaboration based on the first ^1 fb^^ of data taken 
during the 2011 LHC run at 7 TeV. The pMSSM is a 19-dimensional parametrization 
of the MSSM that captures most of its phenomenological features. It encompasses, 
and goes beyond, a broad range of more constrained SUSY models. Performing a 
global Bayesian analysis, we obtain posterior probability densities of parameters, 
masses and derived observables. In contrast to constraints derived for particular 
SUSY breaking schemes, such as the CMSSM, our results provide more generic 
conclusions on how the current data constrain the MSSM. 



1 Introduction 



With the successful operation of the Large Hadron ColUder (LHC) and its detectors in 
2010-11, and with excellent prospects for the future, the LHC is ready to shed light on 
the most pressing open issues in particle physics: the mechanism of electroweak (EW) 
symmetry breaking and the nature of the new physics beyond the Standard Model that 
stabilizes the EW scale (see, for example, Ref. [1]). 

A wealth of theories that extend the Standard Model (SM) have been put forth during 
the past decades. Among these, supersymmetry (SUSY) is arguably the best motivated 
and certainly the most thoroughly studied (see, for example, Refs. [2l[3] for recent re- 
views). Indeed, searches for SUSY rank among the primary experimental objectives of 
the LHC. So far, however, no signal of new physics has been observed at y/s = 7 TeV [4]; 
consequently, the SUSY mass scale has been pushed into the TeV region. 

It is important to note, however, that in the interpretation of their experimental 
results, both the ATLAS and CMS collaborations typically use a very special the- 
oretical model, the so-called Constrained Minimal Supersymmetric Standard Model 
(CMSSM), which is characterized by just four-and-a-half parameters [5l|6]: a univer- 
sal scalar mass mo, gaugino mass mi/2 and trilinear coupling Aq defined at the GUT 
scale Mgut ~ 10^^ GeV, plus tan/3 and sign(/i). The simplifying assumption of uni- 
versality at the GUT scale makes the model very predictive and a convenient showcase 
for SUSY phenomenology. Indeed, it is interesting to present limits within the CMSSM 
because it provides (to some degree) an easy way to show performances, compare limits 
or reaches, etc.. On the other hand, the interpretation of experimental results in the 
(mo,mi/2) plane risks imposing unwarranted constraints on SUSY, as many mass pat- 
terns and signatures that are possible a priori are not covered in the CMSSM. The same 
of course holds true for Simplified Models [7j. 

In this Letter, we therefore present a more general approach, using a 19-dimensional 
parametrization of the MSSM called the phenomenological MSSM (pMSSM) [8]. Us- 
ing results from three independent CMS analyses — the ot hadronic [9j, the same-sign 
dilepton [lOj and the opposite-sign dilepton [11] analyses — we derive constraints on the 
SUSY particles with as few simplifying assumptions as possible. 

The pMSSM parameter space has been thoroughly scanned and studied previously 
in Refs. |12fll5j . It is important to note that [12j "the pMSSM leads to a much broader 
set of predictions for the properties of the SUSY partners as well as for a number of 
experimental observables than those found in any of the conventional SUSY breaking 
scenarios such as mSUGRA [CMSSM]. This set of models can easily lead to atypical 
expectations for SUSY signals at the LHC. " 

The purpose of this Letter is thus to initiate a systematic study that begins with 
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an assessment of what current LHC data tell us, and do not tell us, about the pMSSM. 
We choose to conduct our study using the Bayesian approach [161 |T7] because of its 
conceptual coherence and the direct (intuitive) manner in which probabilistic statements 
are interpreted, namely, as the degree of belief, or plausibility, of a given statement. A 
detailed Bayesian study of the pMSSM was performed in Ref. [13j; however this was 
before LHC data were availableo 

We introduce the pMSSM and its parametrization in Section [21 and outline our 
analysis in Section [3l Our results are presented in Section |H Section [5] contains our 
conclusions. 



2 Parametrization 

The pMSSM, a 19-dimensional realization [S] of the R-parity conserving MSSM with 
parameters defined at the SUSY scale, MsusY = y/^^Wj^H},^ employs only a few plausible 
assumptions motivated by experiment: there are no new CP phases, the sfermion mass 
matrices and trilinear couplings are flavor-diagonal, the first two generations of sfermions 
are degenerate and their trilinear couplings are negligible. In addition, we assume that 
the lightest supersymmetric particle (LSP) is the lightest neutralino, Xi- We thus ar- 
rive at a proxy for the MSSM characterized by 19 real, weak-scale, SUSY Lagrangian 
parameters: 

• the gaugino mass parameters Mi, M2, and M3; 

• the ratio of the Higgs VEVs tan /3 = f 2 /f 1 ; 

• the higgsino mass parameter \x and the pseudo-scalar Higgs mass rriA] 

• 10 sfermion mass parameters m^, where F = Qi, Ui, Di, Li, £^1, Qa, U3, D^, L^, E3 
(imposing rriQ^ = mg^, m^^ = "^Za' ^*^-); and 

• 3 trilinear couplings At, Af, and At , 

in addition to the SM parameters. 

For each pMSSM point, we use Sof tSUSYS .1.6 [19] to compute the SUSY spectrum, 
SuperlsovS.O [20] to compute the low-energy constraints, micrOMEGAs2 . 4 [21] for the 
SUSY mass limits, and HiggsBounds2 . . |22j for the limit on the mas^. Moreover, 
we use SUSYHIT (SDECAYl.Sb, HDECAY3.4) [23J to produce SUSY and Higgs decay 

^In Ref. [18] . the CMSSM gluino-squark mass limits based on 1 fb~^ of LHC data were applied to a 
dark matter global fit in a 9-parameter realization of the MSSM. 

^In evaluating the Higgs mass limit, we apply a Gauss-distributed theoretical uncertainty with a = 
1.5 GeV to the ruh computed by with SoftSUSY, cf. row 8 in Table [T] 
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Table 1: The preLHC experimental results that are the basis of our pMSSM parameter 
scan using Markov Chain Monte Carlo (MCMC) sampling. We re-weight a posteriori 
with the new hmit BR{Bs /u//) < 1.08 x 10"^ at 95% CL However, this has 

hardly any effect. 
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tables, and micrOMEGAs2 .4 [21j to compute the LSP relic density and direct dection cross 
sections. The various codes are interfaced using the SUSY Les Houches Accord [23] . 

3 Analysis 

As noted in the Introduction, the purpose of this Letter is to assess what current data 
tell us, and do not tell us, about the pMSSM. It is convenient to partition these data into 
preLHC and LHC experimental results, which we list in Tables [1] and [21 respectively. We 
use the former to construct a prior 7r{9) on the pMSSM parameter space, which, when 
combined with a likelihood function, L(LHC|0), pertaining to the LHC results, yields the 
posterior density p(0|LHC) ~ L(LHC|0) 7r(0) over the pMSSM parameter space. Here, 
6 denotes the 19 pMSSM parameters Mi, ■ ■ ■ ,At-. We also consider the SM parameters 
nit, rnb{mh) and as{Mz), which are treated as nuisance parameters (see Tabled]). This 
partitioning allows us to assess the impact of the current LHC results on the pMSSM 
parameter space while being consistent with other existing constraints. 
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Table 2: LHC measurements used in the current study. The ot variable is effective 
in suppressing background from hght-quark QCD. SS 2i, and OS 2i denote same-sign 
and opposite-sign dileptons, respectively. The qt [9j, SS [lOj, and OS [11] results were 
published by the CMS Collaboration. 
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The prior vr(^) is constructed as follows. We construct the joint likelihood function 
of the seven independent preLHC measurements D = Di, • • • , D^, of the associated 
observables = lM,---,fJ-7, listed in Table [TJ From Bayes theorem, with a flat prior 
7r(^) = constant □ for each of the seven observables, we obtain the posterior density 
p{f-i\D) = L{D\^) 7r(/i)/p(D) from which we create a random sample of 1.5 x 10'' pMSSM 
parameter points using a standard MCMC technique. During the sampling, we impose 
the constraints on the mass, m/j, of the light neutral Higgs boson (given in row 8 of 
Table [T} and the SUSY mass limits (row 9). Moreover, as we cannot scan over an 
infinite volume, we restrict the sampling to the sub-space |Mj|, |^|,mA, nrip < 3 TeV, 
\At,b,T\ < 7 TeV, and 2 < tan/3 < 6o|l| The MCMC sampling from p{n\D) together with 
the predictions fii = fi{9) induce a distribution over 9 that we take as our prior over 9. 
By construction, the resulting set of pMSSM points are automatically consistent with 
the preLHC experimental constraints listed in Table [TJ 

From the 1.5 x 10^ Markov-chain points we draw a subset of 5 x 10^ points, for each 

^Note that for a Gaussian density, the reference prior is flat (see, for example, Ref. [32] and references 
therein). We will comment on prior dependence in the results section. 

^Evidently, for quantities that are not well bounded by the data within the chosen sub-space, the 
probabilities we calculate will be somewhat sensitive to the choice of sub-space. 
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of which we generate lOK events using PYTHIA6 [33]. (We checked that both the original 
chains and their subsets had converged.) We simulate the response of the CMS detector 
using the publicly available general purpose detector simulation package Delphes [33] . 
Note, that for studies of this scope, a fast, accurate, detector simulation is essential. 
Regarding LHC results, we use the following three published CMS SUSY analyses: 

• the ot hadronic analysis [S], based on 1.1 fb^^, > 2 jets and ax > 0.55, where a-r 
is used to suppress light-flavor QCD, and 8 disjoint bins in Ht, the scalar sum of 
jet transverse momenta; 

• the same-sign (SS) di-lepton analysis [10], based on 0.98 fb~^, with 8 overlapping 
analysis regions of which we use one, Ht > 400 GeV, missing transverse energy 
$T > 120 GeV, and 

• the opposite-sign (OS) di-lepton analysis [llj, based on 0.98 fb~^, with 2 overlap- 
ping analysis regions of which we use one, Ht > 300 GeV and ]^t > 275 GeV. 

We take the observed event counts and background estimates directly from the official 
results of these analyses. For each of the ten results listed in Table [21 we assume a 
Poisson likelihood, 

Poisson ( Nj \sj + bj), ( 1 ) 

with observed count Nj and expected count Sj + bj , where Sj and bj are the expected sig- 
nal and background counts, respectively, for the j'^^ experimental resultj^l Each pMSSM 
point yields predictions for the values Sj,j = 1, • • • , 10. We model the (evidence-based) 
prior for the background parameters bi with a gamma density, 

gammaiK jbjlQj + 1) = e"^^^^ (Kjbjf^ /T{Qj + 1) . (2) 

Here Qj = {Bj/SBj)"^ and Kj = Bj/5Bj, with Bj±6Bj the background estimate in which 
6Bj is taken to be half the width of the confidence intervals listed in Table [21 For each 
pMSSM point, and for each result listed in Table [21 we compute the (marginal) likelihood 
p{Nj\sj) by integrating over the expected background bj. Since, by construction, the 
results are disjoint, the overall LHC likelihood L(LHC|0) is simply the product 

10 

L{L}iC\e) = Y[p{Nj\sj{9)). (3) 

i=i 

The posterior density p(^|LHC) oc L(LHC|^) 7r(0) is approximated by weighting each 
pMSSM point by -L(LHC|^). Finally, we normalize the posterior density over the pMSSM 
sub-space. 

^We use lower-case letters for parameters and upper-case letters for measured quantities. 
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4 Results 



We now present the results of this analysis. Figure [T] shows marginalized 1-dimensional 
(ID) posterior probability density functions of various sparticle and Higgs masses. The 
light blue histograms represent the preLHC probability densities, i.e. taking into account 
only the data listed in Table [TJ Note that the xf and eL/jj-L are bound to be light by the 
Aa^ constraint. Note also that our preLHC distributions differ somewhat from those 
presented in Ref. P3] as we have chosen not to impose any constraint on il/i^. 

The blue, green and red lines show, respectively, the effects of the OS di-lepton, SS 
di-lepton and ot hadronic CMS analyses. The dashed black lines show the final posterior 
densities after inclusion of the results of all three analyses. It is evident that with current 
LHC data-sets, the di-lepton analyses have very little effect on the posterior densities, 
while the ax hadronic analysis pushes the gluino and l*'*/2"'^-generation squark masses 
towards higher values. We also note the slight effect on the Xi LSP mass. The masses of 
other sparticles, including charginos, sleptons and 3^'^-generation squarks, are basically 
unaffected by the current LHC results. This contrasts with the CMSSM case, in which 
all these masses are correlated through their dependence on mi/2 and niQ. Finally, we 
see that the Higgs mass distributions, including that of nih, remain unaffected by current 
SUSY searches. 

The ID distributions of BR{h — )• 57), BR{Bs — )■ Aa^ and the neutralino relic 

density Q.h? are shown in Fig. [2l In Fig. [3] we show the posterior densities of some SUSY 
Lagrangian parameters: M2, ^, At and tan/3. We observe a slight preference for /x > 
with > 0) « 0.53, both pre- and post-LHC startup. This is, however, inconclusive 
(as is the preference for ^ < found in |13)). We confirm the sign correlations between 
Mi and and between At and /i, already demonstrated in Ref. [13]. The corresponding 
plots are available at Ref. [36] . 

It is also instructive to consider correlations between different sparticle masses. Fig- 
ure H] demonstrates the impact of the CMS analyses in the {mu^^rrig) and (m^j^, m^) 
planes. It is interesting to note that the boundaries of the LHC 95% Bayesian credible 
regions (BCRs) approximately match the 95% CL exclusion limits in the CMSSM0 We 
deduce that rrig^q > 1.1 TeV for rrig ~ m^, and rrig > 700 GeV for rriq ^> rrig, is a robust 
conclusion that persists beyond the CMSSM or Simplified ModelsJil 

Our approach moreover allows the study of dependencies between other masses in 

®The way to calculate the BCRs is to some extent a matter of choice. Here we select the region 
containing the highest posterior density values. This is equivalent to choosing the minimal area that 
contains the 68% or 95% of the total volume. 

'^The characteristics of the pMSSM points with < 700 GeV that are not excluded by the current 
SUSY analyses will be the subject of a subsequent study; see also Refs. [151135) in this context. 
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Figure 1: Marginalized ID posterior densities of sparticle and Higgs masses. 




Figure 2: Marginalized ID posterior densities of BR{b — )■ 57), BR{Bs — >■ ^^Ai ), SUSY 
contribution to [g — 2)^, and neutralino relic density VLh?'. 
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Figure 3: Marginalized ID posterior densities for M2, /u, At and tan/3. 
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Figure 4: Marginalized 2D posterior densities of gluino versus squark masses, on the left 
before and on the right after taking the CMS searches into account. The grey and black 
contours enclose the 68% and 95% Bayesian credible regions, respectively. 

a straightforward way, as illustrated in Fig. [S] by means of posterior densities in the 
{m^o,mg) and {m^±,mg) planes. We now see explicitly that bounds on the gluino mass 
are not reflected in chargino and neutralino masses, as would be the case in the CMSSM 
(or actually any scheme with gaugino-mass universality). Moreover, such plots permit 
other interesting observations. In particular, we see how the sensitivity of CMS searches 
to the gluino mass worsens for increasing neutralino or chargino mass. Additional plots 
of ID and 2D distributions are available at Ref. [36j . 

At this point a comment is in order regarding prior dependence. Our preLHC distri- 
butions are of course subject to the same prior dependence that was discussed in Ref. [T3|, 
and persists for quantities that are not much affected by the CMS measurements. How- 
ever, we expect that as the effect of data becomes more influential, the sensitivity to 
the prior diminishes. In fact, for the gluino and l**/2"'^ generation squark masses, the 
likelihood based on LHC data already dominates the prior. For these quantities the 
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Figure 5: Marginalized 2D posterior densities of gluino versus neutralino and of gluino 
versus chargino mass, on the left before and on the right after taking the CMS searches 
into account. The grey and black contours enclose the 68% and 95% Bayesian credible 
regions, respectively. 

posterior distributions are indeed found to be insensitive to the choice of priors. (This 
is also true for some other quantities such as the ii and masses, which are already 
well-constrained by preLHC data.) 

It is also interesting to consider the interplay with other, non-SUSY, searches. Re- 
garding the Higgs sector, in particular the results on H/A — )• rr may have some impact 
on the pMSSM global fit. ATLAS and CMS searches for H/A — )• rr currently exclude 
tan/3 < 10-20 for tua < 250 GeV, and tan/3 < 50-60 for = 450-500 GeV 
While it will be interesting to include this in our global analysis, we note that our 95% 
BCR in the (tan /?, uia) plane, displayed in Fig. [6l starts at tua ~ 500 GeV and shows 
no significant dependence on tan/3. 

Finally, we illustrate in Fig. [7]the interplay with dark matter searches. On the left, we 
show the posterior density in the {VL^h? ,m^o) plane. While matching the cosmologically 
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Figure 6: Marginalized 2D posterior density of tan (3 versus niA- The grey and black 
contours enclose the 68% and 95% Bayesian credible regions, respectively. The preLHC 
distribution looks essentially the same. 

observed value Qh"^ = 0.1123 ± 0.0035 [38j needs a high de gree of fine tuning, a Xi that 
is at least part of the dark matter has a probability of about 60%; see also the right- 
most plot in Fig. [21 In fact there is a slight increase from < 0.123) = 0.53 with 
preLHC data to p{Q-y-h?' < 0.123) = 0.59 when including the CMS analyses, scarcely 
depending on the exact value of the upper bound on Q^h? . The right plot in Fig. [7] 
shows the posterior density of the spin-independent scattering cross section off protons, 
^cr^^(XiP), for the case that the LSP is at least part of the dark matter. Here we imposed 
Q.h? < 0.13 and rescaled the cross section by a factor ^ = r2/i^/0.1123. Note that the 
most credible region is yet to be tested by the direct dark matter searches. 



5 Conclusions 

We presented the first interpretation of the 2011 LHC results based on ~lfb^^ of data 
within the framework of the phenomenological MSSM — a sufficiently generic and well- 
motivated 19-dimensional parameterization of SUSY defined at the SUSY scale. We have 
used three independent LHC SUSY analyses, namely, the CMS ax hadronic, opposite- 
sign dilepton and same-sign di-lepton analyses for this purpose, and expressed our results 
in terms of posterior probability densities. 

Our bounds on gluino and l**/2"'^-generation squark masses match those derived 
in the CMSSM by the experimental collaborations. In addition, we were able to make 
independent statements on the masses and properties of the other SUSY (and Higgs) 
particles, and to show relations between masses that weaken the current bounds. In the 
chargino-versus-gluino-mass plane, for instance, the boundary of the 95% BCR can go 
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Figure 7: Marginalized 2D posterior densities of (left) and rescaled spin-independent 
scattering cross section off protons (right) versus LSP mass. For the latter, only points 
with nh"^ < 0.13 are taken into account. The grey and black contours enclose the 68% 
and 95% Bayesian credible regions, respectively. 

down from rrig ~ 800 GeV to nig w 400 GeV, depending on the xf niass and the rest 
of the spectrum. Our results thus show that current SUSY searches at the LHC provide 
rather limited constraints on supersymmetry in general. With the currently available 
data and searches, we have indeed been able to probe only a small portion of the vast 
pMSSM parameter space, while many regions are still waiting to be explored. Being able 
to work constructively with generic multi-parameter models such as pMSSM will serve 
as a guide to identify the unexplored regions and devise a broader range of dedicated 
searches sensitive to these. 

We have demonstrated in this study that the interpretation of LHC results in terms 
of broad classes of multi-parameter SUSY models is feasible with the currently available 
computational and statistical tools, and that it is indeed possible to make meaningful 
statements on the nature of such models and therefore on supersymmetry, in general. It 
will be interesting to extend our study to include also results from non-SUSY searches. 
Indeed, one of the major advantages of our approach is that it is very well suited for 
global analyses of multiple results from the LHC and elsewhere. 
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